home *** CD-ROM | disk | FTP | other *** search
/ IRIX 6.2 Development Libraries / SGI IRIX 6.2 Development Libraries.iso / dist / complib.idb / usr / share / catman / p_man / cat3 / complib / dlasq1.z / dlasq1
Text File  |  1996-03-14  |  2KB  |  67 lines
  1.  
  2.  
  3.  
  4. DDDDLLLLAAAASSSSQQQQ1111((((3333FFFF))))                                                          DDDDLLLLAAAASSSSQQQQ1111((((3333FFFF))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      DLASQ1 - DLASQ1 computes the singular values of a real N-by-N bidiagonal
  10.      matrix with diagonal D and off-diagonal E
  11.  
  12. SSSSYYYYNNNNOOOOPPPPSSSSIIIISSSS
  13.      SUBROUTINE DLASQ1( N, D, E, WORK, INFO )
  14.  
  15.          INTEGER        INFO, N
  16.  
  17.          DOUBLE         PRECISION D( * ), E( * ), WORK( * )
  18.  
  19. PPPPUUUURRRRPPPPOOOOSSSSEEEE
  20.         DLASQ1 computes the singular values of a real N-by-N bidiagonal
  21.         matrix with diagonal D and off-diagonal E. The singular values are
  22.         computed to high relative accuracy, barring over/underflow or
  23.         denormalization. The algorithm is described in
  24.  
  25.         "Accurate singular values and differential qd algorithms," by
  26.         K. V. Fernando and B. N. Parlett,
  27.         Numer. Math., Vol-67, No. 2, pp. 191-230,1994.
  28.  
  29.         See also
  30.         "Implementation of differential qd algorithms," by
  31.         K. V. Fernando and B. N. Parlett, Technical Report,
  32.         Department of Mathematics, University of California at Berkeley,
  33.         1994 (Under preparation).
  34.  
  35.  
  36. AAAARRRRGGGGUUUUMMMMEEEENNNNTTTTSSSS
  37.      N       (input) INTEGER
  38.              The number of rows and columns in the matrix. N >= 0.
  39.  
  40.      D       (input/output) DOUBLE PRECISION array, dimension (N)
  41.              On entry, D contains the diagonal elements of the bidiagonal
  42.              matrix whose SVD is desired. On normal exit, D contains the
  43.              singular values in decreasing order.
  44.  
  45.      E       (input/output) DOUBLE PRECISION array, dimension (N)
  46.              On entry, elements E(1:N-1) contain the off-diagonal elements of
  47.              the bidiagonal matrix whose SVD is desired.  On exit, E is
  48.              overwritten.
  49.  
  50.      WORK    (workspace) DOUBLE PRECISION array, dimension (2*N)
  51.  
  52.      INFO    (output) INTEGER
  53.              = 0:  successful exit
  54.              < 0:  if INFO = -i, the i-th argument had an illegal value
  55.              > 0:  if INFO = i, the algorithm did not converge;  i specifies
  56.              how many superdiagonals did not converge.
  57.  
  58.  
  59.  
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.